arXiv: 1506.00938vl [astro-ph.HE] 2Jun2015 


Mon. Not. R. Astron. Soc. OOP. fTl fTTR 2015) Printed 3 June 2015 (MN style file v2.2) 


Time series of high resolution spectra of SN 2014J 
observed with the TIGRE telescope 

D. Jack/* M. Mittag/ K.-P. Schroder/’^ J. H. M. M. Schmitt/ A. Hempelmann/ 
J. N. Gonzalez-Perez/ M. A. Trinidad/ G. Rauw^ and J. M. Gabrera Sixto^ 

^ Departamento de Astronomm, Universidad de Guanajuato, A.P. 36000 Guanajuato, GTO, Mexico 

^ Hamburger Sternwarte, University of Hamburg, Gojenbergsweg 112, 21029 Hamburg, Germany 

^ Groupe d’Astrophysique des Hautes Energies, Institut d’Astrophysique et de Geophysique, Universite de Liege, Alice du 6 Aout, 

Bat B5c, 4000 Liege, Belgium 

^ Universidad de Guanajuato, Lascurdin de Retana No 5, G.P. 36000 Guanajuato, GTO, Mexico 


Accepted xxx Received xx; in original form xxx 


ABSTRACT 

We present a time series of high resolution spectra of the Type la supernova 2014J, 
which exploded in the nearby galaxy M82. The spectra were obtained with the HEROS 
echelle spectrograph installed at the 1.2 m TIGRE telescope. We present a series of 33 
spectra with a resolution of i? ~ 20, 000, which covers the important bright phases in 
the evolution of SN 2014J during the period from January 24 to April 1 of 2014. The 
spectral evolution of SN 2014J is derived empirically. The expansion velocities of the 
Si II P-Cygni features were measured and show the expected decreasing behaviour, 
beginning with a high velocity of 14,000 km/s on January 24. The Ca II infrared triplet 
feature shows a high velocity component with expansion velocities of > 20,000 km/s 
during the early evolution apart from the normal component showing similar velocities 
as Si II. Further broad P-Cygni profiles are exhibited by the principal lines of Ca II, 
Mg II and Fe II. The TIGRE SN 2014J spectra also resolve several very sharp Na I D 
doublet absorption components. Our analysis suggests interesting substructures in the 
interstellar medium of the host galaxy M82, as well as in our Milky Way, confirming 
other work on this SN. We were able to identify the interstellar absorption of M82 
in the lines of Ca II H & K at 3933 and 3968 A as well as K I at 7664 and 7698 A. 
Furthermore, we confirm several Diffuse Interstellar Bands, at wavelengths of 6196, 
6283, 6376, 6379 and 6613 A and give their measured equivalent widths. 

Key words: supernovae: individual: SN 2014J - galaxies: individual: M82 - galaxies: 
ISM - ISM: lines and bands. 


1 INTRODUCTION 


Supernovae of Type la (SN la) are of special inter¬ 
est for cosmology since the discovery t hat the expansion 
of the universe i s act ually accelerated jRiess et al.l 1 19981 : 
IPerlmutter et all Il999ll . But des pite being widely used as 
a cal i bratable “standard candle” dPhillipjl 19931: iRiess et all 
1 19961 : IPhillips et al.l Il999l : iGoldhaber et all 200lll this type 
of sup ernovae remains poorly understood. See lParrent et all 
bOlJll for a recent review on SNe la and their properties. 
The physical nature of a Type la supern ova progenitor is 
still under discussion (iLevanon et ahlboiSll as well as the ex¬ 
plosion mechanism. Several differ ent explosion models have 
been suggested and ca l culated llNom otd 1984: Khokhlov 
Il99ll : iPlewa et al.l l2004 [Jordan et al.l l2008l : iPakmor et al 
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l2012l . l2013l : iRosswog et al.l l2009l : iKushnir et al.l l2013l l . A 
well-resolved time series of spectra of SNe la may shed some 
light on the physical properties of the expanding envelope 
and, therefore, might reveal some vital clues about the na¬ 
ture of SN la explosions. 

So far, however, no such detailed time series of SNe la 
spectra exist. The most detailed time series obtained of SN 
2011fe has only a relat ively low spectral resolu tion and some 
gaps in the coverage (IPereira. R. et aObOlSll . The nearby 
and, therefore, quite bright supernova 2014J provided an 
excellent opportunity for observing a SN la very closely and 
obtaining quite high resolution spectra. 

SN 2014J was discovered in the starburst galaxy M82 
llFossev et al.l[20lT during the night of January 21, in an 
early stage of its outbrake. Only two weeks later, around 
February 4, the maximum U-magnitude was reached with 
V = 10.5. By mid-march, SN 2014J had faded below 12.0 
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mag, but it was still well placed in the sky for good-quality 
observations. It has been obser ved by many telescopes and 
in different wavelength r egions IITelesco. C. M. et al.Jl201sl: 
FolewR;_J;_et_aJj 2014: Ashall et al.ll2014l : iMargutti et all 
20141 : Perez-Torres. M. A. et ahT 2014ll . however. none of 


them have obtained observations with a very detailed time 
series of optical spectra with higher resolution. 

As an additional bonus of this bright supernova, we can 
use it simply as a point light source to probe the interstellar 
and perhaps intergalactic medium in the line of sight, mostly 
belonging to the supernova host galaxy and the Milky Way. 
Therefore, high resolution spectra of supernova explosions 
in other galaxies reveal also information about substructures 
in the interstellar medium (ISM) in our galaxy and in the 
respective host galaxy. The lines of interest are those of suit¬ 
able ion species of metals showing resonance lines (e.g. Na I, 
Ca II, Mg II), which represent transitions from the ground 
level, e.g. with a lower level of zero excitation energy. In the 
optical, the sodium D line is most suitable to probe inter¬ 
stellar matter and its dynamics. First high-resolution spec¬ 
tra of this kind were obtained and interpreted on the ex- 
traordinary opportun ity presented by the bright SN 1987A 
dde Boer et al.l [l987tl . Nearly 30 years of advances in de¬ 
tector technology allow such studies to be undertaken with 
much less bright supernovae and more modest equipment. 

Other still unidentified features of interstellar ab¬ 
sorption, which can be observed in SN spectra, are the 
Diffuse Interstellar Bands (DIB). In our Galaxy, such 
DIBs are quite commonly observed in high-resolution 
spectra of massive stars (iHerbid Il995l i. There exist 
also me asurements for the Small and Large Magellanic 
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Cox et al.l I 2 OO 6 I : ICox. N. L. J. et ah. I l2007l : IWeltv et alJ 
. Some spectral observations of extragalactic SNe have 


already succeeded in detecting DIBs (ID ’Odoric o et al.l[l989l : 
ISollerman et al.l[20051 : ICox fc Patadl2008l . l2014l i. For a fur¬ 
ther recent study of DIBs a ga lactic nucleus has been used a s 
the background light source dRitchev fc WallersteirJ[2015li . 
Hence, the nearby SN 2014J presents a further opportunity 
to study DI Bs in high spectral resolution and to expand on 
the work by IWeltv et al.l (l2014h . 

There already exist some work on ISM observations in 
SN 2014J spectra. A detailed analysis of the inte rstellar 
medium of M82 is presented in iRitchev et al.l ( |2015|). They 
use six spectra and determine some abundances. I Weltv et al.l 
ll2014h did a thorough analysis of DIBs using the same set 
of spectra. ICraham. M. L. et al.J d2015l ) present a series of 
very high resolution spectra (i? « 110,000) of SN 2014J. 
They reveal many substructures in the Na I D line. Further¬ 
more, they identify other interstellar absorption lines and 
DIBs. 

In the following section we first introduce the instru¬ 
mentation used, the 1.2 m TIGRE telescope and its HEROS 
spectrograph, then present our observations: a detailed time 
series of 33 high resolution spectra of the Type la supernova 
2014J in M82, covering the period from January 24 to April 
1. We continue with an empirical study of the evolution 
of the prominent Si II line and other spectral features. In 
a further section [3] we take advantage of the high resolu¬ 
tion of the TIGRE/HEROS spectra and study the multi- 
component interstellar absorption in the Na I D absorption 
lines and elsewhere in the SN 2014J spectra. 


2 SPECTROSCOPIC MONITORING OF 

SN 2014J 

2.1 Instrumentation: el TIGRE 

The TIGRE telescope is a fully automated telescope with 
an aperture of 1.2 m, situated near the city of Guanajuato 
in Central Mexico. It is equipped with the HEROS echelle 
spectrograph, which has a resolution of R « 20, 000. Spectra 
are recorded simultaneously in two channels, blue and red, 
covering the large wavelength range from 3800 to 8800 A 
with just a small gap of 130 A around 5800 A. Like opera¬ 
tions, the data reduction pipeline is also fully automatic. 

Originally designed to monitor point-like objects down 
to about 10th magnitude with high-quality spectra, SN 
2014J presented a challenge as much as an opportunity for 
this relatively small telescope. For a more detailed technical 
descript ion of the TIGRE instrumen t ation and its capabili¬ 
ties, see [Schmitt. J. H. M. M. et al, ](l20li). 


2.2 Time series of SN 2014J spectra 

Our monitoring of the supernova 2014J started very shortly 
after its discovery on January 21. We obtained high- 
resolution spectra with a good signal to noise (S/N) of 
around 60 in the red channel in almost every night un¬ 
til March 2. Due to some tecnical problems with the tele¬ 
scope we were only able to take two late time spectra on 
March 31 and April 1. In the blue channel the S/N has 
lower values of around 20, because SN 2014J suffered from 
a significant interstellar absorption and reddening. There¬ 
fore, all spectra presented in Fig. [T] have been dereddened 
with the values E(B — V) — 1.33 and Rv = 1.3 found by 
lAmanullah R. et al.J (l2014l l in their study of the extinction 
law of SN 2014J. Furthermore, since the full resolution is 
not necessary for a study of the broad spectral features of 
supernova ejecta, we binned our spectra to a resolution of 
AA = 10 A for this purpose. Hence, the S/N was further im¬ 
proved and the spectra smoothed, also reducing the visibility 
of some telluric lines inherent to ground based observations. 

Figure [T] shows the complete set of SN 2014J spectra 
observed with the TIGRE telescope. We combined the ob¬ 
servations of the blue and the red channel to present the full 
spectra in the whole wavelength range from 3800 to 8800 A 
except for the small gap of the spectrograph around 5800 A. 
The dates at the right hand side mark the observation date 
in 2014 in Universal Time (UT). Because of the strong inter¬ 
stellar reddening of SN 2014J, the later spectra in the blue 
channel have a poor S/N and one cannot clearly distinguish 
features any more. 

The very first spectrum was observed on January 24 UT 
and, therefore, still about 10 days before maximum light. 
This spectrum already shows the prominent P-Gygni profile 
of the Si H feature at around 6300 A, which usually identifies 
a SN as a Type la. A high expansion velocity feature of the 
Ga H triplet is also present in the early spectra of SN 2014J 
at a wavelength of around 8000 A. We were able to observe a 
spectrum almost every night until March 2 UT. In this time, 
31 spectra were obtained. Then, due to a simple technical 
problem with the telescope mount, which caused some out- 
time for repair work, we could only obtain two more late¬ 
time spectra of SN 2014J, on March 31 and April 1 UT. 







































































Time series of SN 201 fJ spectra with TIGRE 3 



4000 


5000 


6000 7000 

wavelength [A] 


8000 


24 

Jan 

25 

Jan 

26 

Jan 

27 

Jan 

28 

Jan 

30 

Jan 

31 

Jan 

01 

Feb 

02 

Feb 

03 

Feb 

04 

Feb 

05 

Feb 

08 

Feb 

12 

Feb 

13 

Feb 

14 

Feb 

15 

Feb 

16 

Feb 

17 

Feb 

18 

Feb 

19 

Feb 

20 

Feb 

21 

Feb 

23 

Feb 

25 

Feb 

26 

Feb 

27 

Feb 

28 

Feb 

01 

Mar 

02 

Mar 


Figure 1. Spectral time series of SN 2014J observations obtained with the TIGRE telescope in both the red and blue channel of the 
HEROS spectrograph. The dates are given for the observations in early 2014. 
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wavelength [A] 

Figure 2. Evolution of selected SN 2014J spectra obtained in the 
blue channel of the HEROS spectrograph. 


A SN la spectrum at early times is supposed to be quite 
flat and should not show many clear features. This is clearly 
consistent with the appearance of our early-time spectra of 
SN 2014J. For clarity, the latest spectrum shown in Fig.[T]is 
the one observed on March 2, which corresponds to about 1 
month after maximum light. The small ’’lines” in some parts 
of the spectra stem from a few, more noticeable telluric lines. 

The here presented time series reveals the changing 
shape, on an almost daily schedule with some features. The 
general trend, the P-Cygnl extrema moving towards the red¬ 
der part of the spectrum, is consistently observed in the Si II 
line and some other features. This means that the observed 
expansion velocities change during the evolution. The reason 
for that is the expansion of the envelope, which decreases the 
mass densities and, therefore, decreases all opacities. Hence, 
with ongoing expansion, ever deeper and slower layers of the 
expanding envelope become visible in the spectra. In other 
words, the quasi-photosphere moves inwards into less fast 
expanding shell material. 

Another strong feature in SN la spectra is caused by the 
Ca II IR triplet at around 8500 A which already appears 
during the early phases at around maximum light. It can 
be observed throughout the whole evolution until the latest 
spectrum shown in Fig. [T] At early times it shows a high 
velocity component. 

At later phases an Fe II emission feature appears about 
where the prominent Si II line was observed earlier. Another 
Fe II feature then appears at around 5000 A. Indeed, it is 
true of evolved type la SNe to show more spectral features in 
general. Eventually, these go into emission, as the late, opti¬ 
cally thin phase of the expanding envelope has been reached. 
At the same time, one sees the inner layers, which are rich 
of iron peak elements, causing the above-mentioned rise of 
Fe II features. In the following section, we present a more 
detailed picture of the evolution of selected spectral features. 

2.3 Detailed spectral evolution of important 
featnres 

In Fig. m we compare a selection of only four spectra from 
the blue channel of the HEROS spectrograph, representative 


of four different epochs in the evolution of SN 2014J. Due 
to the strong interstellar reddening of the light of SN 2014J, 
the later spectra have only poor S/N. 

The first spectrum from above in Fig. [2] was obtained 
on January 25. It shows the typical features of a SN la spec¬ 
trum during maximum light. We can see the typical S H 
line profile in form of a W between wavelengths of, roughly, 
5000 and 5500 A. At shorter wavelengths a clear Fe H fea¬ 
ture can be seen at around 4500 to 5000 A. Furthermore, 
around wavelengths of about 4000 to 4500 A a clear fea¬ 
ture of Mg H is present in the spectrum of SN 2014J. Below 
4000 A one can see the drop in flux towards the Ca H H 
& K feature, although it is not fully covered in the shown 
spectrum and the S/N is already quite low in that part of 
the spectrum. We should stress that spectra of SNe la are a 
blend of millions of lines, so that one cannot always assign 
one specific element to each observed feature. 

The spectrum of February 4, which corresponds to max¬ 
imum light, shows similar spectral features, when compared 
to the spectrum of January 25. However, the line profiles 
have already changed a bit. Additionally, this spectrum 
shows consistently lower expansion velocities, in all of its 
features. The further two spectra shown in Fig. [2] of Febru¬ 
ary 18 and the one of March 1, suffer already from poorer 
S/N. Nevertheless, these spectra show the expected emission 
features of Fe H. 

2.3.1 Evolution of the Si II line profile 

The broad Si H P-Cygni line profile is shaped by the ejecta 
of the supernova, its dynamics and its opacity evolution, in 
a most representative way, since this is the strongest feature 
observed in the SN la spectra at 6355 A during maximum 
light. 

A time series of the wavelength range of this Si H fea¬ 
ture is shown in Fig. [3] A set of 10 spectra was selected, 
which best represent the evolution and its chronology here. 
All these spectra are shown in their full spectral resolution. 
They are normalized in flux, and their wavelength scales 
were transformed onto a barycentric velocity scale (with re¬ 
spect to 6355 A). 

In the first spectra, the Si H feature exhibits a typical P- 
Cygni profile, which is expected for a rapidly expanding, op¬ 
tically still thick envelope. Hence, the emission is still much 
less prominent than the absorption part of the profile. Fur¬ 
thermore, it is clearly seen from the individual line profiles, 
how fundamentally the absorption part changes. At the be¬ 
ginning, this is just one consistent broad absorption, as of a 
typical P Cygni profile. Its minimum is located at an expan¬ 
sion velocity of around ~ —14, 000 km/s. This characteristic 
expansion velocity then shifts slightly to « —11,000 km/s 
until the day of February 20. By February 26, however, the 
Si H absorption has changed its nature significantly, because 
a small emission feature appears inside of the Si H P-Cygni 
absorption feature. By April 1, the Si H feature has come 
close to disappearance. 

During the later phase an Fe H emission feature arises 
in the wavelength range of the Si H feature. During the on¬ 
going expansion the envelope becomes thinner and allows 
that deeper parts of it shape the spectra. In this way, the 
iron peak elements from the inner envelope of the expand¬ 
ing SN la envelope become visible. By the day of February 
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Figure 3. The Si II line of 10 observations are shown. These 10 spectra are shown the chronology evolution of the time series Si II line 
at 6355 A. 


26 a small emission feature at around « —7000 km/s has 
appeared in the absorption trough of the P-Cygni profile 
of the Si II feature. The emission feature to the right at 
« +9000 km/s is an Fe II feature, which becomes stronger 
until April 1. By that time the envelope has become optically 
thin. 


2.4 Expansion velocities of Si II and Ca II 

Since we obtained a well resolved time series of SN 2014J 
spectra, we were able to study the expanding envelope in 
some detail. For that, we measured the evolution of the ex¬ 
pansion velocities of the Si II feature at around 6300 A and 
the Ca II IR triplet feature at around 8500 A. 
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Figure 4. Measured expansion velocities of the Si II feature at 

6000 A. 



Figure 6. Full resolution capture of the complex IS absorption 
in the sodium D-line doublet, from the superposition of the 32 
best spectra. 
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Figure 5. Measured expansion velocities of the Ca II feature at 
8500 A. During the early phase it shows a high velocity compo¬ 
nent. 


Figured] shows the change in the expansion velocity of 
the prominent Si II feature. Starting from a high expansion 
velocity of 14,000 km/s on January 24, the expansion ve- 
locites decrease while the envelope is expanding. Hence, the 
deeper parts are revealed in the spectra, which have slower 
expansion velocities. When the Si II is about to disappear, 
the expansion velocity is about 12,000 km/s. The photo¬ 
sphere has then passed the layers with abundant silicon. 

In Fig. (5] we show the measured expansion veloci¬ 
ties of the Ca II feature at around 8500 A. Thanks to 
our early observations, we could study the evolution of 
a high velocity component of the Ca II feature. Its ex¬ 
pansion velocity decreases from 24,000 km/s in our first 
spectrum down to about 20,000 km/s in the spectrum 
6 days lat er. In their study o f near infrared spectra of 
SN 2014J, iMarion. G. H. et al.J ll2015ll found a similar be¬ 
haviour for the Ca II IR triplet expansion velocities. They 
measured values for the high expansion velocity feature be¬ 
tween 26,000 km/s and 20,000 km/s. Like in our study, they 
did not detect a clear high velocity component of the Si II 
feature. All these observations are consistent with the study 


of high ve locity components in ot her type la supernova per¬ 
formed bv ISilverman et ahl ll2015ll . Observing a high velocity 
component of Ca II is much more common than one of Si II. 

The main component of the Ca II feature starts with an 
expansion velocity of roughly 14,500 km/s and decreases to 
a value of 12,000 km. As expected in SNe la, the observed 
expansion velocities of the Ca II feature are similar to those 
measured for the Si II feature. 


3 THE INTERSTELLAR ABSORPTION 

COMPONENTS OE THE SODIUM D-LINES 

As can be seen in Fig. [6l the D-line doublet of Na I at 
5889.95 A and 5895.92 A shows a complex, multi-component 
absorption line structure due to distinctly different contri¬ 
butions. The telluric absorption lines were not removed in 
the presented spectra. We here study and discuss this very 
valuable information and offer an interpretation in terms of 
different origins. Within the noise of individual spectra, the 
D-line absorption features appeared invariable over the two 
months or so we succeeded with the spectroscopic monitor¬ 
ing. To minimize noise, we superimposed the 32 best SN 
spectra in the vicinity of the Na D doublet (see Fig. ID), af¬ 
ter a barycentric correction was applied to the laboratory 
wavelength scale of each spectrum. 

There are essentially four absorption components or 
groups of even narrower absorption features, as have also 
been see n by other groups in e ven h i gher spectral resolu ¬ 
tion (seelGraham. M L. et al.J ll2015^ : IWeltv et ^ ll2014h : 
iRitchev et al.1 I ^OlSh i: (a) is a small absorption (the right¬ 
most feature in Fig. EJ at a radial velocity shift of -1-245 km/s 
(barycentric), which is clearly seen in the longer-wavelength 
Na I D-line but is blended (see below) in the left one. (b) 
is a very broad and saturated absorption trough, promi¬ 
nent in both D-lines. Its distinct, asymmetric profile in both 
D-lines suggests an interesting velocity and density struc¬ 
ture, which reaches from a sharp edge at -1-50 km/s to be¬ 
yond -1-200 km/s, close to component (a). Component (c) is 
a sharp, medium-strong component well visible in both D- 
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lines at —5 km/s, while (d) is a 4 times weaker (by equivalent 
width) component, equally sharp and distinct, at —60 km/s. 
It is well visible only for the shorter-wavelength D-line, while 
its longer-wavelength counterpart (d’) falls into component 
(a) of the shorter-wavelength D-line. 

The uncertainties of the equivalent widths, as measured 
in the integrated spectrum are about 0.015 A for the weak 
lines (components a, c and d, see below) and 0.05 to 0.1 
A for the broad component b. Individual spectra yielded 
uncertainties about 6 times as large. Consequently, the weak 
components were well resolved but much affected by noise. 
Any hidden temporal changes, if any, would have had to 
be smaller than about 25% of the absorption of the weak 
components. 

3.1 Components (c) and (d): distinct ISM clouds 

For component (c) the obvious absorber candidate is galactic 
ISM. The line of sight to SN 2014J probes high galactic lat¬ 
itudes away from the galactic centre, in galactic coordinates 
toward s I = 141.4° and b = 40.6°. Around this galactic lon¬ 
gitude, [Moha^^eF^i] (l2004l l observed radial velocities of cold 
ISM clouds between —12 and +10 km/s (see their Fig. 2), 
with respect to the local standard of rest (LSR). A small 
correction of +6 km/s applies to translate these veliocities 
into the heliocentric system quoted above. Still, component 
(c) fits well in here. 

Despite its modest appearance, component (c) suffers 
from some degree of saturation, because the two versions 
of it do not show equivalent widths in a ratio of 1:2 (as of 
the respective f-values) but differ by only a factor of 1.25 
(0.20 and 0.16 A, ±0.01 A, approximately). Furthermore, 
their line widths are very close to the instrumental profile, 
which suggests a very small broadening velocity much under 
10 km/s. Therefore, saturation can set in already with small 
equivalent widths. 

The analysis of pure absorption lines was formerly ap¬ 
plied by us to special problem s in t he field of stellar at¬ 
mospheres. See lSchroeder et al.l lll994l ~) and references given 
there for a brief summary of the simple theory. Given well- 
known transition probabilities (for th e Na I D-line doubl et 
we use the / and g values given bv IWiese et al.l lll969lD . 
the equivalent width of a weak, unsaturated line is a linear 
function of only the column number density N of the ab¬ 
sorber. However, once saturation sets in, at least two lines 
of the same multiplet and with much different transition 
probability / are required to solve for both, N and the line- 
broaden ing velocity. We here u se the same simple computer 
code as ISchroeder et aP lll994ll . to calculate the profiles of 
pure absorption lines and their equivalent widths for any 
such parameter pair. 

In the case of the two Na I D (c) components, we find 
a best-matching solution for A^Na = 3.3(±0.7) x 10^^ cm~^ 
and Udisp = 3.5(±0.2) km/s. T his small velocity dispersion 
matches the lower end of what IWilson. C. D. et al..l (l201ll l 
(see their Table 1) find for whole galaxy disks, and presum¬ 
ably we indeed see just a single cloud in our own Milky 
Way. We should note, however, that a possibly smaller ve¬ 
locity dispersion naturally allows for a larger column de nsity 
due to saturation effects. In fact. I Ritchey et aP (l2015l l find 
9.5 X 10^^ cm“^ at 2.0 km/s for this component. The sodium 
column density only gives us a lower limit for neutral hydro¬ 


gen of Nm = 2.1 (±0.5) X 10^® cm“^, when using a relative 
sodium abundance of 10“®’®, since IS sodium is much more 
prone to ionization (much lower ionization potential) than 
hydrogen. The tru e value may be almos t two orders of mag¬ 
nitude higher fsee lRitchev et all (l2015lj '). 

Apparently, the distinct, small and sharp component 
(d) is of a similar nature (single cloud) as component (c). 
Even with the same small broadening velocity, it would not 
suffer much from saturation, as the equivalent width is only 
0.05 A. Accordingly, for its column density we obtain A^Na ~ 
3 X 10^^ cm“^. At —60 km/s barycentric radial velocity, it 
may appear not very likely to be of a galactic origin, but the 
Leiden/Argentine/Bonn (LA B) H 21cm Survey do es show 
galactic IS hydrogen here (see iRitchev et all ll2015h L 

3.2 The signatures of cool ISM in the M82 disk 

Now, can the remaining D-line absorption components (a) 
and (b) be attributed to the M82-internal ISM? And is the 
observed D-line absorption consistent with the IS extinc¬ 
tion in the line of sight towards SN 2014J, which mainly 
arises in the disk of M82 (as the line of sight stays clear 
of most of the galactic IS absorption)? The maximum mag¬ 
nitude of SN 2014J of V = 10.5 mag suggests Ay ~ 1.9 
mag (with M — m = 27.8 mag at 3.6 Mpc distance and 
My,max = —19.2 mag for a typical SN type la). The large 
reddening (see above) seems to suggest even a larger IS ab¬ 
sorption. _ _ 

Using the relation given by iGiiver fc Ozell ll2009ll . i.e. 
AH/(cm“^) = 2.2 X 10^^Av/(mag), a neutral hydrogen col¬ 
umn density of A^h of about 4 x 10^^ cm“^ should be ex¬ 
pected. The respective sodium column density would be lim¬ 
ited to A^Na < 6 X 10^® cm“^. Due to ionization, the true 
value could be nearly two orders of magnitude lower. 

The (here clearly not given) case of an optically thin 
absorption translates the equivalent width of all D-line ab¬ 
sorption, nearly 4 A, into a lower limit for A+ja of only 
« 2 X 10^® cm“^. Hence, this is indeed fully consistent with 
an IS absorption of about 2 mag (or more). Nevertheless, 
from a prohle analys is of a much better spectral resolution, 
IRitchev et al.l (l20lki hnd only 1.8 x 10®'* cm~^. 

M82 as a whole h as a radial velocity of ±203 km/s 
llChvnoweth et al]l2008f) . With respect to any Na D absorp¬ 
tion components, however, only ISM, which lies in front 
of SN 2014J, becomes visible. SN 2014J is located in the 
western disk-half, 58” away from the galaxy nucleus. Here, 
the ISM rotates towards us wit h (in the rest-frame of M8 2) 
Wrad ~ —100 to —150 km/s fsee lMavva fc Carrasco! ll2009lH . 
For this reason, the D-line absorption should appear between 
(barycentric) radial velocities of ±50 km/s to ±100 km/s. 
Remarkably, the broad component (b) features a sharp blue 
edge at ±50 km/s, but then reaches out with a red wing to 
component (a) at ±245 km/s (barycentric). That material 
is clearly not participating in the general disk rotation but 
(in the rest-frame of M82) moves backwards, towards the 
SN 2014J. 

An answer to such diverse dy namics may be give n by 
the hot superbubble discovered bv iNielsen et al.l ll2014h . Its 
X-ray emission appears to coincide with the site of SN 2014J. 
The bubble has been inflated to about 200 pc in diameter 
by the action (hot winds and SNe) of some very massive 
stars formed there perhaps 50 Myrs ago. We suggest that 




















































8 D. Jack et al. 



Figure 7. Detection of absorption by ISM for the Ca II 3933 and Figure 8. Detection of absorption by ISM for the K I and K I 

Ca II 3968 lines in the high resolution spectra of SN 2014J. 7698 lines in the high resolution spectra of SN 2014J. The three 

narrow lines at the left hand side are telluric lines. 


the fully saturated part of component (b) between +50 and 
about +140 km/s (barycentric) traces the swept-up, cold IS 
gas in front of that superbubble. In the rest-frame of M82, 
it moves with about —150 km/s to —60 km/s. Part of this 
gas is pushed towards us by the still expanding superbub¬ 
ble, in addition to be rotating towards us with the disk of 
M82. In particular, this would explain the sharp blue edge 
of component (b): this is the cool material closest to the 
bubble-shockfront approaching us. 

The backward moving ISM found in the red wing of 
component (b) and in component (a) could be explained, if 
SN 2014J was not located inside the superbubble, but in¬ 
stead a bit behind it, placing the superbubble at +150 km/s 
(barycentric), that is, in co-ration with the disk, and have 
its shockfront still expanding at about 100 km/s. In this 
scenario, SN 2014J would also probe some ISM pushed 
out away from us, towards it and so moving against disk 
rotation. The small but distinct component (a) at about 
+40 km/s in the rest-frame of M82 would mark this cold 
ISM closest to the far-side bubble-shockfront. However, the 
hi ghly complex substructure of the IS absorption resolved 
bv iGraham. M. L. et al.,] ll2015ll leaves, as intriguing as our 
proposal is, alternative explanations, e.g. by a number of 
individual structures in the long line of sight. And compo- 
nent (a) may als o ma tch the western H I s treamer found by 
IYuu et all J 19931 ) , see iRitchev et al.l ll2015l) . 

3.2.1 Other IS absorption signatures: K I, Ca 11 and 
DIBs 

Interstellar absorption, mostly of component (b) belonging 
to the galaxy M82, is visible in several other lines as well, 
foremost in the Ca II doublet. Figure [7] shows two absorp¬ 
tion features corresponding to Ca II H and K at wavelengths 
of 3933.66 A and 3968.5 A. Since these lines are at the lower 
end of the blue channel, the spectra at that wavelength 
range have low S/N. Obviously, one cannot do a very de¬ 
tailed study. However, both absorption features are broad, 
ranging between +50 and +200 km/s, much like component 
(b) of the D-lines of the Na I doublet. We were also able to 
detected two ISM absorption feature of K I at 7698.96 A and 




Figure 9. Two features of DIBs observed in the high resolution 
spectra of SN 2014J. The left hand side shows the feature at 
6196 A, and on the right hand side the DIB 6613 A. 


at 7664.9 A as shown in Fig. [S] Note that telluric absorp¬ 
tion lines have not been removed from the observed spec¬ 
tra. These feature also show the broad absorption at around 
+ 100 km/s corresponding to the broad IS absorption fea¬ 
ture of the sodium D-line doublet. This broad feature is due 
to absorption of ISM in M82. We were also able to detected 
CH+ at a wavelength of 4232 A. 

The high resolution of our spectra allows us to check 
even for unidentified interstellar absorption features, the so- 
called Diffuse Interstellar Bands (DIB). So far, the responsi¬ 
ble elements or molecules have not been identified, which is 
a motivation for further high-quality observations of DIBs. 
We were able to detect a few of them in our best high reso¬ 
lution spectra of the supernova SN 2014J. 

Since the DIB features are present in all our spectra, we 
again added them up to obtain just one final spectrum with 
further improved S/N, As shown in Fig. [9] we find two DIBs 
at wavelengths of 6196 A and 6613 A. The DIB wavelengths 
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radial velocity [km/s] radial velocity [km/s] 

Figure 10. Three more features of DIBs observed in the high 
resolution spectra of SN 2014J. We found DIB features at 6283 A, 
as shown on the left hand side, and on the right hand side two 
features at 6376 and 6379 A. 


Table 1. Measured equivalent widths of the ISM and DIB fea¬ 
tures observed in the spectra of SN 2014J. 


Feature Equivalent Width [A] 


DIB 6283 

0.27 

± 

0.04 

DIB 6376 

0.065 

± 

0.013 

DIB 6379 

0.037 

± 

0.007 

DIB 6613 

0.2 

± 

0.01 

K I 7664 

0.35 

± 

0.05 

K I 7698 

0.28 

± 

0.03 


have been taken from the study of iHobbs. L. M. et al.J 
ll2009h . Again, the plot indicates velocities in the range of 
roughly +50 km/s to +200 km/s, as found in the broad (b) 
absorption feature of the Na D doublet. 

In this combined spectrum of SN 2014J, we also found 
several other DIBs. Figure [TOl shows the DIB of 6283 A and 
the double feature of DIBs at 6376 A and 6379 A, in the 
same range of +50 k m/s to +200 km/s . All these feature 
were also observed bv lWeltv et al.l (l2014l l and several other 
authors from high resolution spectra of SN 2014J. We were 
unable to detect the 5780 and 5797 A DIBs since these fea¬ 
tures unfortunately lie exactly in the wavelength gap of the 
HEROS spectrograph. 

In Table [U we present the measured equivalent widths 
of the ISM and DIB features observed in our spectra of 
SN 2014J. The spectra in the wavelength range of the Ca II 
H & K absorption lines are too noisy to able to obtain 
reasonable equivalent widths. Our measured values for the 
eq uivalent width s of th e DIBs are in the same range a s foun d 
by IWeltv et al.l ll2014l ~) and lOraham. M. L. et al.J ll2015h . 
However, for the DI B 6283 A, the stro ngest of the observed 
DIBs in our spectra. I Weltv et al.l (l20I4h found a much higher 
value for its equivalent width, while Graham. M. L. et al..l 
(I 2 OI 5 II found a value comparable to our. The reason might 
be the broad wing s of this line as seen in the spectra of 
IWeltv et al.l (|2014l l. In our spectra, the telluric lines make 
it difficult to measure equivalent widths including the broad 


wings (see Fig. nsi). The IS absorption in both lines of K I 
are stronger than the DIBs. 


4 CONCLUSIONS AND DISCUSSION 

We presented a time series of 33 high resolution spectra of 
the Type la supernova 2014J in M82 observed with the TI¬ 
GRE telescope. We obtained these spectra with the HEROS 
echelle spectrograph, which covers the wavelength range 
from 3800 A to 8800 A, with only a small gap at around 
5800 A. The long monitoring period of about two months, 
well including maximum light of SN 2014J, in combination 
with the high quality of the spectra, enables us to empiri¬ 
cally describe the spectral evolution of this supernova la in 
detail. By its spectral features and their behaviour SN 2014J 
is very typical for a Type la SN. 

During maximum light, red-channel SN 2014J spectra 
show broad P-Cygni profiles caused by the strongest lines 
of Si II and Ca II. In the blue channel we observed the typ¬ 
ical features of S H, Mg II, Fe II and Ca H. We studied 
the evolution of the prominent P-Cygni profile features of 
Si H at around 6300 A in some detail. It shows a decrease in 
the expansion velocity in time as it is expected from rapidly 
expanding SN envelopes. The first spectrum from January 
24 shows a minimum in the absorption trough at an expan¬ 
sion velocity of ~ —14,000 km/s. The expansion velocity 
decreases, and on February 20 it has dropped to a value of 
« -12,000 km/s. 

The Ca H infrared triplet feature at about 8500 A shows 
a high expansion velocity component during the first days of 
observation with expansion velocities of 24,000 km/s mea¬ 
sured on January 24 decreasing to 20,000 km/s right be¬ 
fore this component disappears. The main expansion veloc¬ 
ity component shows a similar behaviour as the Si II feature. 
The expansion velocities decreases from an inicial value of 
14,500 km/s down to 11,0000 km/s during the later phases 
of the expansion of the envelope. 

Due to the significant amount of interstellar absorption 
and reddening of SN 2014J, the later blue-channel spectra 
suffer from a quite low S/N. In the well-exposed red channel, 
however, those later spectra document the transition from 
the optically thick to the thin phase of the expanding en¬ 
velope very well. Now, emission features are dominated by 
Fe II lines from deeper layers. At the same time, in the red 
channel, the Si H feature disappears, while a Fe II feature 
arises in the respective wavelength region. 

The reason for the rise of these Fe II lines lies in the lay¬ 
ered structure of the SN debris, combined with a decreasing 
opacity: the quasi-photosphere moves inwards and passes 
from the outer layers of intermediate mass elements like sil¬ 
icon and sulfur to deeper layers with iron peak elements, 
stemming from the inner envelope of the exploded progen¬ 
itor White Dwarf. In this way our dense time series allows 
us to basically scan the abundance properties of the entire 
SN envelope, especially its abundances structure. Achieving 
this requires, of course, a huge amount of modeling work. 
Every spectrum of the series needs to be analyzed and/or 
modeled with a suitable (ID or even 3D, non-LTE, dynamic) 
atmosphere code . In a future work, using th e PHOENIX code 
llHauschildtII I 992 I : iHauschildt fc Baronlll999lf , we will do this 
work for the whole time-series of spectra, in oder to obtain 
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a detailed and quantitative abundance structure of the en¬ 
velope of SN 2014J as a case study for a SN of Type la. 
Since the different explosion models predict slightly different 
abundances and also a differently layered abundance struc¬ 
ture, this detailed abundance determination might exclude 
or favor one or some of the suggested explosion mechanism 
for a supernova of Type la. Some of the here presented spec¬ 
tra have already been used to identify the feature, which 
causes the secondary maximum in the 7-band light curve of 
Type la supernovae djack et al.ll2015l ). 

Furthermore, the high resolution of 7? « 20, 000 of the 
obtained spectra allows for a study of the sharp ISM ab¬ 
sorption features in the SN 2014J spectra. The Na I D dou¬ 
blet shows different substructures in the ISM of our galaxy, 
as well as of the host galaxy M82. There is a small fea¬ 
ture of ISM of M82 at a high velocity of about -1-245 km/s 
(barycentric). A very broad absorption feature at velocities 
of -1-50 km/s to +200 km/s is apparently caused by ISM 
in M82, perhaps even related to the expanding environment 
around a superbubble just in front of the SN. We also found 
two sharp absorption features at velocities of about —5 km/s 
and —60 km/s, both are of a galactic origin. We are also able 
to identify the same broad ISM absorption of M82 between 
+50 km/s and +200 km/s in the principal lines of Ca II H 
& K at 3934 and 3968 A and K I at 7664 and 7699 A. 

Thanks to the high quality of the spectra, we also iden¬ 
tified a few DIBs in the spectra of SN 2014J, confirming 
other groups results. We found signatures of these features 
at wavelengths of 6196, 6283, 6376, 6379 and 6613 A and 
measured the equivalent widths of these feat ures. Hence, 
we ca n fully confirm the observations made bv lWeltv et al.l 
(l2014h in their analysis of high resolution SN 2014J spectra. 
The radial velocities of the observed DIBs again are fully 
consistent with the broad Na D absorption feature of M82 
in the range of velocities of +50 km/s to +200 km/s. 
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